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The dynamics of flexible polymers in dilute solutions is studied taking into account the hydrody- 
namic memory, as a consequence of fluid inertia. As distinct from the Rouse-Zimm (RZ) theory, the 
Boussinesq friction force acts on the monomers (beads) instead of the Stokes force, and the motion 
of the solvent is governed by the nonstationary Navier-Stokes equations. The obtained generalized 
RZ equation is solved approximately. It is shown that the time correlation functions describing 
the polymer motion essentially differ from those in the RZ model. The mean-square displacement 
(MSD) of the polymer coil is at short times ~ t^ (instead of ~ t). At long times the MSD contains 
additional (to the Einstein term) contributions, the leading of which is ~ y/i. The relaxation of the 
internal normal modes of the polymer differs from the traditional exponential decay. It is displayed 
in the long-time tails of their correlation functions, the longest-lived being ~ t~^^'^ in the Rouse 
limit and t^^''^ in the Zimm case, when the hydrodynamic interaction is strong. It is discussed that 
the found peculiarities, in particular an effectively slower diffusion of the polymer coil, should be 
observable in dynamic scattering experiments. 



PACS numbers: 36.20.Ey, 82.35.Lr, 83.80.Rs, 05.40.Jc 



I. INTRODUCTION 



Most of the theoretical investigations on the dynamic 
properties of flexible macromolecules performed so far are 
based on the Rouse-Zimm (RZ) model HH H 0. In 
this theory the polymer molecule is modeled as a chain 
of beads under Brownian motion. The bead participates 
in the interactions with solvent molecules. The solvent 
contributes a frictional force against the motion of a 
monomer and a random force to take into account the 
random collisions exerted on monomers. In the Rouse 
model, the solvent is considered nonmoving. Within the 
Zimm theory, the motion of each monomer also affects 
other monomers, via the flow which it induces in the am- 
bient medium. It has been shown that the Zimm model 
predicts the correct dynamical behavior for dilute poly- 
mer solutions in ^-conditions. The Rouse model is appli- 
cable for good solvents, where the corrections due to the 
hydrodynamic interactions and excluded volume effects 
cancel each other to a large extent, or in the situations 
when the surrounding polymers screen out the hydrody- 
namic interactions |3, Q] . Both models assume Gaussian 
equilibrium distribution of the beads, which is in most 
cases the only description of a polymer that can be ma- 
nipulated analytically. The models hold for such poly- 
mer properties, which involve length scales that are large 
compared to monomer sizes. Although the RZ model is 
considered a universal theory well describing the long- 
time behavior of the polymer macromolecules, there is 
still a number of unresolved problems in the understand- 
ing of the polymer dynamics in solution |^|^. So, the 

dependence of the first cumulant in the intermediate 
scattering vector regime is well confirmed for synthetic 



polymers and DNA, however, the experimental values 
are smaller than the theoretical predictions. Systematic 
deviations from the theoretical behavior at large scatter- 
ing vectors q have been found for various polymers using 
quasi-elastic neutron scattering. At small q the diffusion 
coefficient calculated from the continuous RZ model also 
deviates from the experimental values for both the natu- 
ral and synthetic macromolecules. For a review of these 
and other difficulties of the RZ theory we refer also to 
the earlier work ; it can be seen that the situation has 
changed little during a decade. 

In this work we propose a generalization of the RZ 
theory that could provide a better understanding of the 
dynamical behavior of polymers in dilute solutions. The 
bead and fluid inertia during the motion of the polymer 
in the solvent have been taken into account. Also the 
hydrodynamic interaction is considered by solving the 
nonstationary Navier-Stokes equations. A possible im- 
portance of such a development (in connection with the 
Rouse model) was already discussed in Ref. 8. As dis- 
tinct from the traditional approach to the polymer dy- 
namics 0, S 01 J in our theory the resistance force 
on the moving bead will not be assumed the Stokes one, 
which is proportional to the velocity of the bead. We 
model this force by the Boussinesq-Basset history force 
0,0,^ which, at a given time t, is determined by the 
state of the bead motion in all the preceding moments of 
time. We thus have a possibility to obtain solutions valid 
for almost arbitrary t (for incompressible fluids, however, 
t >> b/c, where 6 is the monomer size and c the sound 
velocity). The discussed generalization means that the 
effects of hydrodynamic memory are taken into account. 
Such effects are extensively studied in the physics of sim- 
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pie liquids and in the theory of Brownian motion (see e.g. 
Ref. jl^l) with very important consequences. In partic- 
ular, the memory effects are revealed in such interest- 
ing peculiarities like the famous " long-time tails" of the 
molecular velocity autocorrelation function (VAF), first 
discovered by means of computer experiments [T^l IT3 | . 
The concept of the Brownian motion lies in the basis of 
the RZ theory of polymer dynamics. Thus it is natural 
to expect that the memory effects are important to poly- 
mers as well. An attempt to show it was done in our 
recent paper [3| where the Rouse model has been con- 
sidered taking into account the hydrodynamic memory 
but with no hydrodynamic interaction. Here we present 
a more general theory that, as limiting cases, includes 
both the Rouse and Zimm models. It will be shown that 
the inclusion of the hydrodynamic memory into the clas- 
sical RZ model leads to an essentially different behavior 
of the time correlation functions describing the polymer 
motion in solution. It will be demonstrated by the ap- 
pearance of long-time tails of these functions that reflect 
a strong persistence of the correlation with the initial 
state of the polymer. The relaxation of the mean square 
displacement of the whole polymer, as well as the decay 
of the polymer internal modes, are slower than in the 
original model. We also show how the tails are displayed 
in the long-time behavior of the Van Hove function used 
in the interpretation of quasi-elastic scattering of light 
and neutrons. It is discussed that the found new fea- 
tures in the polymer dynamics should be observable in 
the scattering experiments. Moreover, we believe that 
the obtained results could contribute to the solution of 
some of the existing problems between the theory and 
experiment. 



II. THE DYNAMICS OF POLYMERS WITH 
HYDRODYNAMIC MEMORY 

Within the RZ model the motion of the nth polymer 
segment (the bead) of a flexible polymer coil consisting 
of N beads is described by the equation 



M 



(1) 



Here, Izf is the position vector of the bead, M is its mass, 

7ch 
„ is the force from the neighboring beads along the 

chain, 7„ is the random force due to the motion of the 

— ^fr 

molecules of solvent, and / „ is the friction force on the 
bead during its motion in the solvent. In the RZ model 
the latter force was 



7'^ 

J n 



d'X'r, 



dt 



(2) 



with if (lzf„) being the velocity of the fluid in the place 
of the nthe bead due to the motion of other beads. The 
friction coefhcient for a spherical particle of radius b is 
f = Qmjb, where 77 is the solvent viscosity. However, 



this expression holds only for the steady-state flow. In 
the general case the resistance on the body moving in a 
liquid depends on the whole history of the motion, i.e. 
on velocities and accelerations in the preceding moments 
of time. For incompressible fluids we use the Boussinesq 
force H [13 instead of Eq. (2). This means that we 
replace, in the Fourier transformation (FT) with respect 
to the time, the friction coefficient ^ with the frequency 
dependent quantity 



(3) 



where x = y^—iujp/ri, (Rex > 0) and p is the density of 
the solvent. The Boussinesq force follows naturally from 
the usual hydrodynamics as the solution of linearized 
Navier-Stokes equations 0, 0|. Additionally to the 
Stokes force it contains terms which, if the fluid den- 
sity is comparable to the density of the bead, cannot be 
neglected for the nonstationary motion when they are of 
the same order as the inertial term in Eq. (1). Equations 
(1-3) have to be solved together with the hydrodynamic 
equations for the velocity of the solvent. 



dlf 



-Vp + r/Alf + Vlf = 0. 



(4) 



Here p is the pressure. The quantity ip is an external 
force per unit volume jij. 



(5) 



Equations (4) are solved using the FT in coordinates and 
time. The solution can be, for any of the components 
a — x^y, ov z written in the form 

«o (^) = j'^^'T. - V'^ (^') ' (6) 



/3 

with the FT of the Oseen tensor 

Hap (^) = A^ap + Br^rpr^'^, 



(7) 



A = (8^7?r)-^ {e-y - y [{l - e"^) y-']"}, 
B = (STTTyr)"^ {e-y + 3y [(l - e"'^) y-^] "}. (8) 

Here the prime denotes the differentiation with respect 
to y ^ rx- Substituting tp'^ from the FT of Eq. (5) to 
from (6), and the obtained result into the FT of equation 
of motion (1), we get a generalization of the RZ equation, 
which in the continuum approximation reads 

- = ^ in) + /- (n) + Mu^x-^ (n)] (9) 



N 



dmH: 



aPnm 



SksT d^x'^ 
a? dm? 



+ (m) + Muj'^xf^ (m) 
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where a is the mean square distance between neighboring 
beads along the chain. It has been used that the force 
between the beads can be obtained from the effective 
potential u = (3fcBr/2a^) ^^^2(^11 ~ ~^ n-xf' which 
follows from the equilibrium (Gaussian) distribution of 
the beads 0, B- Due to the dependence of the Oseen 
tensor on the difference ~r' nm = '^(n) — lS'{m), Eq. (9) 
is nonlinear and thus hardly solvable analytically. We use 
the common approximation of preaveraging of the tensor 
over the equilibrium distribution P(r„m) = (2TTa'^\n — 



m|/3)-3/2 exp[-3r2^/(2a2|„ _ m\ 



(10) 



h'^ {n — m) 



y^BTT^jn — m\r]a 



with z = xci (1?^ — '7i|/6)^^^. In the case without mem- 
ory [4] the function h at large |n — m| behaves as ^ 
\n — •m\~^l'^\ now the effective interaction between the 
beads disappears more rapidly, ^ |n — m|~^/^. Since 
Eq. (9) now contains only the diagonal terms, it can 
be solved using the FT in the variable n, Izf^ (n) — 
'y^o + 2X]p>i ifp cos{7Tnp/N), where the boundary con- 
ditions at the ends of the chain have been taken into ac- 
count dl^/dn = at n = 0, TV. The inverse FT then 

yields the following equation for the Fourier components 

— >ij 



1 — y/iTze'^ erfc(z) 



i'pit) = {yap{o)yap{t)): 



duj- 



*Re5t^ 



(14) 

in agreement with the FDT. The generalized susceptibil- 
ity is Upiuj) = [{2- Spo)N]-^[-iujE'^ - Mu;"^ + Kp]-^ , and 
the generalized forces corresponding to the coordinates 
Up^ are Nf^^. Using the Kramers-Kronig relation [T^ . 
the same initial value for the function ipp ai t = as 
in the RZ theory is obtained: ipp{0) = kBTap{0) = 
kBTi2NKp)-\p>0. 

Equation (14) gives the solution of the model. Know- 
ing ilip{t), other correlation functions of interest can 
be found, e.g. the VAF, $p(t) = {vap{Q)vap{t)) = 
—(P'il}p{t) I dt^ ^ or the mean square displacement (MSD), 
{Ay^{t)) = 2['0p(O) - i^pit)]. The previous RZ results in 
the absence of memory are obtained by putting lo = 
in Sp, Eq. (11), and M = 0. The mode yo describes 
the motion of the center of inertia of the coil |3, 
In the RZ case we get "00(0) — "00 (t) = Dct. The 
diffusion coefficient Dc — ksT {h^Q + l/N^^ contains 

the Zimm {Dc = SksT (sVeir^riaj ^) and Rouse 

{Dc = ksT/N^) hmits. The internal modes {p > 0) 
relaxed exponentially, V'p(i) = {kBT/2NKp)exp{—t/Tp), 
with the relaxation times Tp = [Kp (l + 2NhppS,) 



where h^p — 



l!% - fp - Mu? 



(11) 



A. The Rouse limit with memory 



where St; 



r 1 



^p^NKv\ 



and 



?,Ti^p^kBTl{Na)\ p = 0,1,2,... The matrix is de- 
fined by the integral 



1 

iV2 



N 



N 



an / amn (n — m) cos cos . 

Jo ^ ' N N 

(12) 

In obtaining Eq. (11) the fact that the nondiagonal el- 
ements of the matrix are small in comparison with the 
diagonal ones and can be in the first approximation ne- 
glected has been already taken into account; the sub- 
stantiation of this is the same as in Refs. [3,4]. Equa- 
tion (11) can be investigated as it is usually done in the 
theory of Brownian motion. One can use the fluctuation- 
dissipation theorem (FDT) or the properties of the 
random forces J^j. The forces acting on different beads 
n and m are uncorrelated, thus their correlator is ~ Snin- 
In going to the continuum approximation the Kronneker 
symbol has to be replaced by the (5-function, S{n — to), 
so in the FT we have 



I^UJ ruj' 
J pa J qo 



kBTReE' 



(2 - Spo) 7:N 



^5ap5pq5{uj +U}'). (13) 



Equation (11) then yields the following expression for the 
time correlation function of the Fourier components yap, 



Now let us consider the case with memory, the limit 
of Rouse |15| . It assumes that the hydrodynamic inter- 
action contribution to in Eq. (11) is negligible for all 
ujs. The subsequent equations then change only by the 
substitution « . The corresponding integral in Eq. 
(14) is encountered in the theory of Brownian motion 
of one particle: in the case when p = Q, the particle is 
free, and if p > 0, it is in a harmonic field with the force 
constant Kp. Such problems were solved in a number of 
investigations beginning from the work |T^ . Adopting 
the known solutions, see e.g. Ref. [12], for the MSD of 
the coil we have the following asymptotic expression: 



{Ayl {t)) = 2Dct 



1 - 



< 9 V MsJ t ' 
(15) 

where t >> Tf, = b^p/rj and Ms is the mass of the sol- 
vent displaced by one bead. For small times (Ai/q {t)) w 
kBTt^/N {M + Ms/2). (The physically correct value 
kBTt^ /N M can be obtained only if the compressibility is 
taken into account ^ISiilfi]-) It is well seen from Eq. (15) 
and confirmed by numerical calculations of Eq. (14) how 
slowly the previous result is approached: the second (Vt) 
term is less than 1 per cent of the first only for t > lO'^rfc. 
For long chains when the terms ~ M can be neglected, 
at i = lOrh the MSD constitutes only about 3/4 of the 
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Einstein limit. Note that for the motion of the coil as 
a whole the exact analytical solution exists which differs 
from the known solution for one Brownian particle only 
by a factor 1/iV TE,'?^. For example, the VAF of the 
center of mass of the Rouse coil is 



$o(0) 
Ai — Aq 



J2 (-l)'+'A,e(^?*)erfc(A,Vt), (16) 



where $o(0) = ksT / {M + Ms/2)N and A^ arc the (com- 
plex) roots of the equation A^ + y^A/r + 1/t = with 
T = (M + M,/2)/e 

The long-time asymptote of the function ipp>o{t) de- 
scribing the relaxation of the Rouse modes is 



1 



Mt) ^ 



5/2 



(17) 



of \J —iio. The first correction to the classical result is 
determined by the coefficient at the term (— iw). The 
first nonvanishing correction to the susceptibility ap(w) 
is given by the term (— iw)^/^. The coefficient at this 
term is determined by the expansion coefficients at the 
terms (— iw) and {—iiSf'l'^ in the expansion of h^^. For 
the correlation function tpp^oit), using Eq. (14), we thus 
find 



V'p(o) 



457r-^ 



16 TR 



1 3/2 



(20) 



where t » tr and Tp = {Na^^^^{r]/kBT){3TTp^y^^^ is 
the Zimm relaxation time 0, Q . 



C. End-to-end vector and the dynamic structure 
factor 



B. The Zimm case 

For the Zimm model, when the hydrodynamic inter- 
action is strong for all frequencies that significantly con- 
tribute to the studied correlation functions, we have in 

"ppi 



Eq. (11) w [{2- Spo)Nh'^p] \ The Oseen matrix 
(12) can be calculated with any degree of precision, e.g. 
for p = we have the exact result 



4:Z 



00 



1 



z \ 



erfcz — 1 



(18) 



where z = (iVxa/6)^/^. Then in the ^-representation 
the correlation function can be obtained using the stan- 
dard methods Here we show the main terms of the 
asymptotic expansion of the time correlation functions 
that can be found also by using the expansion of the sus- 
ceptibility ap(ti') in small (— zw)-'^/^. In the case of the 
diffusion of the coil as a whole we obtain the expression 
similar to Eq. (15), 



V'o(O) - Mt) - Dc 



t - — {TRt) 



1/2 



(19) 



The characteristic time in this equation is tr — R^p/rj. 
The form of Eq. (19) exactly corresponds to that fa- 
miliar in the theory of Brownian motion of rigid par- 
ticles of radius R. Here R is the hydrodynamic ra- 
dius of the coil 0, 01 , determined from the relation 
Dc = ksT / [QirRrj) where Dq is the Zimm diffusion co- 
efficient. In the standard model of Gaussian chains the 
radius R is given also by the Kirkwood relation '5| R~^ = 

iV-2EliELi,™^„(''nm> « SVV^/a, where w 
is the distance between the beads and >> 1. 

Consider now the internal modes of the polymer, p > 0. 
The components of the matrix (12) are calculated as in 
Ref. [4] and can be expressed through special functions, 
so that their expansion is known to any desired power 



Having the solutions (17) and (20), the evolution 
of the end-to-end vector of the chain can be investi- 
gated. From the relation lt{t) = lt{t,N) — lzf(t,0) = 
-4Ep=i,3,...l^pW one finds 



$(t) = (i?(t)i?(0))-48^i(0) 



p=l,3, 



p2^p(0) 



(21) 



for both the models. For the Rouse model at long times 



4 1/2 / o 2 

TT T\T, , , s I on Tl 



(22) 



and in the Zimm case the long-time asymptote reads 
$(t) w -3V2^Vi(0)(rfl + l.85Ti)T^j/^t-^/^. (23) 

(In Eqs. (22) and (23) Tp^i is the relaxation time for the 
corresponding model.) 

Finally, we give the result for the intermediate scat- 
tering function G{kt) — A^"^ X]nm('2^P{^ ^ [^"(^) ~ 
'afm(O)]}) that is used in the description of the dynamic 
light or neutron scattering from a polymer coil pj] ( k 
is the wave- vector change at the scattering). Acting in 
a similar way as in Ref. [3], G{k ,t) can be for large t 
approximated by the expression 



X exp 



G{k,t) 



36 



1 



Arexp{-fc2 [^o(o)- Vo(t)]} 



37r6 



E 

P=2,4,. 



(24) 



valid for kR << 1 (in the opposite case, the function 
G{k, t) becomes very small at long times). Equation (24) 
is equally applicable for both considered models, if the 
corresponding Eqs. (15) and (17) for the Rouse case, or 
(19) and (20) for the Zimm model are substituted here. 
One can see that the contribution of the internal modes 
is small and thus hardly detectable against the diffusion 
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term given by the first exponential. However, our pre- 
dictions concerning the diffusion of the coil as a whole, 

G{k,t)^Nc^p{-k^[MO)-Mt)]}, (25) 
could be directly measured in the scattering experiments. 

III. CONCLUSION 

We conclude that in the generalized RZ model, when 
the memory of the viscous solvent is taken into account, 
the relaxation of the correlation functions describing the 
polymer motion essentially differs from the original the- 
ory. The MSD is at short times ~ (instead of t). At 
long times it contains additional (to the Einstein term) 
contributions, the leading of which is ^ \/t. The inter- 
nal normal modes of the polymer motion now do not 
relax exponentially. It is well displayed in long-time 
tails of their time correlation functions. The longest- 
lived contribution to the correlation function of the bead 
displacement is ^ t^"^/'^ in the Rouse limit and t"'"^/^ in 
the Zimm case, when the hydrodynamic interaction is 
strong. It would be interesting to investigate the found 
peculiarities using computer simulation methods and ex- 
perimentally, e.g. by the dynamic light or neutron scat- 
tering. Simple estimations show that the relaxation of 
the internal modes, ^pyo{t), although qualitatively dif- 
ferent from the previous RZ model, for real polymers 
only slightly differs from the traditional exponential law 
~ exp(— </Tp), except at high frequencies (it is because 
the relaxation time for the internal modes, Tp, is much 
larger than the characteristic times Tb and tr). However, 
the differences from the original model, at least for the 
Zimm case, should be experimentally accessible. Due 
to the long-range character of the hydrodynamic field, 
the characteristic time of the Zimm model, tr, is de- 
termined by the size of the whole coil. For a typical 
radius around 100 nm, the density and viscosity of wa- 
ter at room conditions, one gets tr about 10 nanosec- 
onds {tr/tp^i « 5 X 10""* sec). Taking into account 
the possibilities of current experiments (e.g. in Ref. 
the dynamic structure factor of polymers in solution was 
studied using quasi-elastic light scattering in a time win- 
dow beginning from 12.5 ns, and even shorter times are 
accessible by the neutron spin-echo technique H^), and 



the fact that the function '0o(O) — ipoit) approaches the 
Einstein limit Dct very slowly as t increases, the non- 
diffusive (" ballistic" ) motion of the polymer coil should 
be readily observable by the dynamic light and neutron 
scattering. In fact, similar experiments were successfully 
carried out on single Brownian particles. For example, 
using the diffusive wave spectroscopy the ballistic motion 
of polystyrene spheres with the radius b = 0.206/im in 
aqueous solution (with the characteristic time tr about 
0.04/isec) was observed The size of such particles 

corresponds to the hydrodynamic radius of the DNA coil 
of a molecular weight 6 x 10^ g/mol (with the diffusion 
coefficient L»c « 1-3 x 10"® cmVs) The nondiffu- 

sive motion of even smaller (with a radius < 100 nm) 
particles was observed in the experiments psf . It has 
been found in these and other works that the apparent 
diffusion coefficient of Brownian particles is smaller than 
that following from the Einstein theory. These observa- 
tions are very similar to the situation described by us for 
the polymer coils. The tails in the MSD lead to a slower 
decay of the dynamic structure factor (25). This corre- 
sponds to diffusion with an effectively smaller diffusion 
coefficient (as well as the first cumulant) than predicted 
by the previous theory. This is the long-standing un- 
resolved "puzzle" between the theory and experiments. 
We believe that the proposed theory could help to solve 
this problem in the description of the dynamic scattering 
experiments on polymers 0, 0, Q a-nd thus to contribute 
to a deeper understanding of the dynamical properties of 
polymers. 
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